The Information Geometry 
of the Spherical Model 

W. Janke 
Institut fur Theoretische Physik 
Universitat Leipzig 
Augustusplatz 10/11 
D-04109 Leipzig, Germany 

D.A. Johnston 
Department of Mathematics 
Heriot-Watt University 
Riccarton 
Edinburgh, EH14 4AS, Scotland 
and 
R. Kenna 

School of Mathematical and Information Sciences 
Coventry University 
Coventry, CV1 5FB, England 

February 1, 2008 

Abstract 

Motivated by previous observations that geometrizing statistical mechanics offers 
an interesting alternative to more standard approaches, we have recently calculated 
the curvature (the fundamental object in this approach) of the information geometry 
metric for the Ising model on an ensemble of planar random graphs. The standard 
critical exponents for this model are a = —1, (5 = 1/2, 7 = 2 and we found that the 
scalar curvature, 1Z, behaves as e -2 , where e = (3 C — [5 is the distance from criticality. 
This contrasts with the naively expected 1Z ~ e~ 3 and the apparent discrepancy 
was traced back to the effect of a negative a on the scaling of 1Z. 

Oddly, the set of standard critical exponents is shared with the 3D spherical 
model. In this paper we calculate the scaling behaviour of 1Z for the 3D spherical 
model, again finding that 1Z ~ e~ 2 , coinciding with the scaling behaviour of the 
Ising model on planar random graphs. We also discuss briefly the scaling of 1Z in 
higher dimensions, where mean-field behaviour sets in. 



1 The Information Geometry of Spin Models 



The idea of endowing the space of parameters with a metric and geometrical structure has 
been borrowed from parametric statistics [|lj and employed to some effect in statistical 
mechanics P, §, £|, |, ||, [/], §|. The approach seems to be particularly fruitful for a spin 
model in field where the parameters are (3, the inverse temperature, and h, the external 
field. In this case the (Fisher-Rao) metric is simply given by 

Gij = didjf , (1) 

where / is the reduced free energy per site and d. t = (d/d/3, d/dh). 
For such a metric the scalar curvature may be calculated as 



R = 1 



2G 2 



djf d p d h f dlf 
d}j dp h f dpdlf 
d}d h f dpdlf %f 



(2) 



where G = det(Gjj) is the determinant of the metric itself. 

The work in 0, [| f|, [5| has made it clear that, as one might expect, the scalar curvature 
plays a central role in any attempt to look at statistical mechanics from a geometrical 
perspective. In particular, for all the models that have been considered so far, the curva- 
ture diverges only at a phase transition point for physical ranges of the parameter values. 
We assume a standard scaling form for the free energy per site 

f(e,h) = \- 1 f(e\ a %h\^) , (3) 

where e = j3 c — f3 (with f3 c marking the critical point) and a e , ah are the scaling dimensions 
for the energy and spin operators. 

The spherical model has a one-sided critical point || and we are interested in the 
high-temperature domain (where e > 0). There, standard scaling assumptions allow us 
to write (Q) as 

f(e,h) = e l/a ^+(he- ah/a ^) , (4) 

where we have introduced the scaling function ip + . Consideration of the behaviour of the 
components of 1Z in equ. (0) then leads to the scaling 

n~e a ~ 2 , (5) 

for the curvature itself. 

Equation ([5]) shows that the scalar curvature not only characterises phase transition 
points as divergences, as do the more standard statistical mechanical quantities such as 
the specific heat C and susceptibility %i ft displays a scaling behaviour that allows the 
extraction of the critical exponent a. Indeed, with hyperscaling, equ. @ can be recast as 

, (6) 

where £ is the correlation length and d is the dimensionality of the system. So the 
curvature is proportional to the correlation volume, an intuitively reasonable result on 
dimensional grounds. 



At first sight, our calculation of the scaling behaviour of TZ in [T(| for the Ising model 
on planar random graphs puts a spanner in the works. There, a = — 1 but TZ ~ e~ 2 rather 
than the expected TZ ~ e~ 3 . However, returning to the detailed scaling of the individual 
terms in equ. (0) showed that a negative a affected some of the scaling of the components 
and the end result could be traced back to these modifications. In principle this modified 
behaviour should apply to any model with negative a, so an interesting test would be to 
calculate TZ, or at least its scaling limit, for another model with this property. 

The spherical model provides just such a test case. It was solved (in field) in the classic 
Berlin and Kac paper |9| and the critical exponents in 3D found to be a = — 1, (3 = 
1/2, 7 = 2 which are identical to those of the Ising model on 2D planar random graphs 
This is remarkable, because there are no obvious physical similarities between the 
two models. An additional motivation for any such calculation is the general paucity of 
spin models that have been solved in field. Indeed, TZ for spin models has been obtained 
explicitly only for the ID Ising ||, Bethe lattice Ising || and ID Potts |12j models and 



the scaling form calculated for the Ising model on planar random graphs |T0|. This forms 
a rather small sample on which to formulate hypotheses about its general behaviour and 
its scaling properties. In the sequel, we extend this short list incrementally by looking 
at the scaling behaviour of TZ in the spherical model. We concentrate on the 3D case, 
but also discuss other dimensions, including the mean-field like behaviour which sets in 
at d = 4. 



2 The Spherical Model 

Berlin and Kac |{J introduced the spherical model (and the Gaussian model) in an attempt 
to understand how generic some of the features of Onsager's solution of the 2D 
Ising model are for ferromagnetic spin models, particularly for other dimensions. In the 
spherical model, the ±1 condition on the value of the Ising spins is relaxed, whilst retaining 
a global constraint on the total spin magnitude. With Si denoting the value of a spin at 
a site « of a hypercubic lattice, the partition function is M 



Z = dsi... ds N exp (3 ^ SiSj + Si \ 5 lJ2 s "i ~ N ) 



(7) 



where N is the total number of sites. This can be evaluated by exponentiating the 
constraint and using steepest descent, resulting in the following expression for the reduced 
free energy per site in the thermodynamic limit, iV — > oo: 

/ .l log (l) +/3 ,_i sW + _^_ , (8) 

where 

1 r 27T ( d \ 

9( z ) = / dux... du d log \z - V cos(cu fc ) . (9) 



The saddle point value of z which appears in the expression for the free energy in equ. (§) 
is determined from 



The solution reveals no transition for d = 1 and 2, and a transition with the exponents 
a = — 1, (3 = 1/2,7 = 2 for d = 3. For ci > 4, on the other hand, mean-field behaviour 
with a = 0, /3 = 1/2,7 = 1 se ^ s m (modified by multiplicative logarithmic corrections in 
the d = 4 case |TJJ). For the following formulae we confine our attention to the d = 3 
case. 

It is useful to note that, with h — 0, equ. ( [TDD gives 

- = ^ (11) 

dp g"{z) ' 1 ] 

and hence 

^ = _4^) 

e^ 2 ~ [g"{z)f • 1 J 

The critical point is given by z = d = 3 and fe = |, and the behaviour of g(z) in this 
region is determined by differentiating equ. (|9|) twice and then expanding for the small 
Uk values which give the dominant contribution. One finds 

s " w ~-^- 3) ~ 1/2 ■ (13) 



A further differentiation gives 



and an integration yields 



^-^-3)-^ , (14) 



^) = 7^-3) 1/2 + s'(3) , (15) 



where g'(3) = (18 + 12v^2 - lOVZ - 7^6) [K(2y/3 + - 2^ - 3)] 2 « 0.505 462 019 .. . 
is the exactly known massless 3D lattice propagator at the origin, with K{k 2 ) denoting 
the standard elliptic integral. This latter expression can be combined with equ. ([[(]) with 
h = to give 

(^-3)~87r 2 (/? c -/?) 2 ~e 2 , (16) 

in which j3 c = g'{3)/2 « 0.252 731 009 . . .. Equations (|TJ) and (|TJ may then be substi- 
tuted in equs. (|lj) and (|i~2"D to give the scaling of dz/df3 and d 2 z/d(3 2 , 

£s# = l^- 4 ^- 3 ) I/2 >=° ■ 



which we shall employ below in the calculation of the scalar curvature. 



3 The Scalar Curvature 



We recapitulate the general considerations of |10] before going on to discuss in detail the 
results for the spherical model. In (|4]), we define A = l/a e and C = —ah/a € , which, in 
terms of the standard exponents, gives A = 2 — a and A + C = (3. \{ A> 2 (which is the 
case for a < 0) the specific heat, rather than diverging, will be a constant at the critical 
point, which we denote by A(A — 1)0(0). With this in mind, substitution of the scaling 
function into equ. (0) gives 

n- 1 



2G 2 

A(A — 1)0(0) e A+2C ^l(Q) 

-A(A-l)(A-2)e A - 3 ^ + (0) -(A + 2C)e A+2C - 1 ^(0) 

-(A + 2C)e A+2C - 1 ^l(0) 

(18) 

In equ. ( |T8"D terms with an odd number of h derivatives have been set to zero. These 
do not appear because of the h — > —h symmetry of the free energy. The scaling for the 
determinant of the metric G is given by 

G = A(A-l)e A+2C m^l(0) , (19) 

so the leading term in 1Z (for A > 2) is 



(A + 2C) 2 
2A(A- 1)0(0)' 



or, translating back to the standard critical exponents, 



7? = 1 f - 2 (oW 

2(2-a)(l-a)0(O) ' 1 } 

In summary, if a < the expected scaling of 1Z is 1Z ~ e~ 2 rather than the TZ ~ e a ~ 2 , 
which is seen for positive a. 

We now move on to examine the scaling of the various terms contributing to TZ in 
equ. (0) for the spherical model itself. As we have remarked, the h — > —h symmetry in the 
free energy per site, /, of the spherical model means that any terms with an odd number 
of h derivatives will automatically be zero when h = 0, hence fp^ = fpph = fhhh = 0. 
This leaves the non-zero terms (again when h = 0) 

hp = ^ + 



fhh 



d(3 2(3 2 
1 

2/3(z-3) 
d 2 z 1 



d(3 2 (3 3 ' 

f 1 9z _ 1 (22] 

Jhhp 2(3{z-3) 2 d(3 2(3 2 {z-3) ' 1 ; 



which, using equs. ([16]) and (|17]) have the following behaviour in the scaling region: 



2fl 



f 1 e- 2 

Jhh ~ 16vr 2 /3 c (/3 c - /3) 2 ~ 



Z/3/3/3 ~ 16vr 



2 1 



Aft/3 ~ 8 ^(3 c ((3 c -(3r 16^(^-^)2 ~ e 3 • (23) 

Comparing with the individual terms in equ. ( TT8T ) we see that the expected general scaling 
of each term (for a < 0) does indeed apply and that overall we have, as in equ. (f2ll), 

K ~ e~ 2 . (24) 

We thus see that calculating the scaling of 1Z for the 3D spherical model for which a = — 1 
gives results in accordance with expectations from general scaling arguments which take 
into account the negative a, similarly to the Ising model on planar random graphs. 

In d > 4 dimensions the exponents of the spherical model attain their mean-field 
values, a = 0, (3 — 1/2,7 = 1- I n this case g'{z) ~ z — d + g'(d) near criticality and the 
resulting scaling of the components of the expression for 1Z is 



n 



2G 2 



const e 1 
const e -2 
e- 2 



(25) 



where the determinant of the metric scales as G = e _1 . Thus 1Z itself scales as e~ 2 , in the 
mean-field case. This is in agreement with the expected scaling of e a ~ 2 when a > since, 
in that case [|J, 

k. 1 



2G 2 

A(A — l)e A ~~ 2 ip + (0) e A+2C ^+(0) 

-A(A-l)(A-2)e A ~ 3 ij + (0) ~(A + 2C)e A+2C ~ 1 ip'^(0) 

-(A + 2C)e A+2C ~V+(0) 

(26) 

where the scaling of the metric determinant is 

G = A(A - l)e 2A+2C - 2 ^j + (0)^l(0) . (27) 



The difference in scaling between the a > and the a < cases originates in the 
term which contributes to both TZ and G in equ. ( p6]) and equ. (p7|). In the negative a 
case, the contribution of this term comes from the regular part of the free energy, and is 
non-diverging. 

Note that for a = 0, A = 2, and although the displayed scaling term A(A — I) (A — 
2)e A ~ z ip + {Qi) in equ. (^6|) would be expected to vanish, a regular term could still contribute. 



This is indeed what is seen in the explicit calculation for the spherical model. Two further 
caveats arise. Firstly, for the upper critical dimension d = 4 itself, there are multiplicative 
logarithmic corrections |14| which require a careful handling. Secondly, for d > 4, where 
the scaling is described by mean-field theory, the hyperscaling relation fails. So one should 
be quite cautious in comparing the explicit calculation in terms of e with expressions 
involving the correlation length £. 



4 Conclusions 

We have calculated the scaling behaviour of the scalar curvature, 7Z, of the information 
geometry metric for the 3D spherical model, finding, as for the Ising model on planar 
random graphs, that TZ ~ e~ 2 . Careful considerations of the scaling of the various elements 
which contribute to TZ show that this is in accordance with expectations. We also briefly 
discussed the scaling behaviour of TZ for d > 4, corresponding to mean-field exponents. 

The calculation reported here provides another example of a statistical mechanical 
model in which the curvature of the information geometry metric diverges at the critical 
point, with a clearly quantifiable scaling behaviour. Once again it is curious that the 
relevant exponent is shared with the Ising model on planar random graphs. All critical 
exponents now coincide in the two models though there is no obvious physical relation 
between them. It seems that the suggestion in the original Berlin and Kac paper || 
that the 3D spherical model might provide a hint to the behaviour of the 3D Ising model 
applies rather to the Ising model on 2D planar random graphs. 
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